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INTRODUCTION 

It  has  been  shownw  that  it  is  possible  to  invert  the  magnetic 
field  first-derivative  equations  for  a  point-magnetic  dipole  field, 
therebv  obtaining  the  bearing  and  scaled  magnetic  moments  (quantities 
related  to  the  magnetic  moments)  of  the  field  source  from  a  knowledge  of 
the  first  field-derivatives  at  a  single  point  only.  This  inversion  of 
the  first-derivative  dipole  equations  can  be  done  in  closed  form  and  is 
greatly  facilitated  by  carrying  out  the  analysis  in  the  principle  axis 
frame  of  the  first  field-derivative  tensor.  Unfortunately  there  seem 
to  be  no  analogous  special  reference  frames  which  simplifv  the  analvsis 
of  the  higher  derivative  equations  for  the  corresponding  higher  multi¬ 
pole  sources.  However,  by  using  a  spherical  tensor  representation, 
wherein  onlv  the  independent  quadrupole  moments  and  field  derivatives 
appear,  it  is  possible  to  formulate  the  inversion  of  the  magnetic  field 
second-derivative  equations  for  a  point  quadrupole  source  in  terms  of 
two  simultaneous  transcendental  equations  which  appear  to  be  amenable 
to  numerical  solution  on  a  digital  computer.  These  equations  have  been 
solved  numerically  for  several  special  cases.  In  each  such  special  case 
it  was  found  that  there  were  not  more  than  six  solutions  and  that  these 
solutions  exhibited  a  certain  degree  of  symmetry  about  the  origin.  If 
it  is  assumed  that  this  is  true,  in  general,  for  the  quadrupole  inver¬ 
sion  problem  it  would  indicate  that  there  are  three  solutions  in  each 
half  space,  making  a  total  of  five  ghost  solutions  and  one  physical 
solution.  This  is  not  substantially  worse  than  the  dipole  inversion 
problem  if  the  half-space  in  which  the  source  resides  is  known  a  priori. 

The  spherical  tensor  procedure  has  also  been  shown  to  be  applicable 
to  a  monopole  source  (for  heuristic  purposes  only)  and  a  dipole  source. 

A  recipe  for  inverting  the  dipole  problem  by  the  spherical  tensor  pro¬ 
cedure  is  given  in  the  appendix.  The  spherical  tensor  procedure  is  also 


''Naval  Ship  Research  and  Development  Laboratory  Report  3493,  Dipole 
Tracking  with  a  Gradiometer ,  by  W.  M.  Wynn,  January  1972, 
Unclassified. 

<'S/ Naval  Coastal  Systems  Laboratory  Report  135-72,  Inversion  of  the 
Magnetic  Field  Gradient  Equations  for  a  Magnetic  Dipole  Field,  by 
C.  P.  Frahm,  November  1972,  Unclassified. 


believed  to  be  applicable  to  the  magnetic  field  Jlth-derivative  equations 
for  a  point  2^-pole  source  although  a  complete  proof  has  not  yet  been 
established.  In  this  report  the  general  scheme  of  the  procedure 
believed  to  be  applicable  to  an  arbitrary  2^-pole  source  will  be  pre¬ 
sented  followed  by  the  details  of  the  inversion  for  a  quadrupole  source. 


MULTIPOLE  FIELD  DERIVATIVE  EQUATIONS 


In  a  static  situation  the  magnetic  field  H  produced  by  a  magnetic 
source  can  be  derived  from  a  magnetic  scalar  potential  $  by  the  familiar 
relation 

3  *■  -V$  (1) 

Furthermore,  if  one  is  willing  to  introduce  the  fiction  of  magnetic 
poles,  the  potential  can  be  conveniently  expressed  in  the  form 


where  p(r')  is  the  magnetic  pole  density  of  the  source  at  point  r'  and 
r  is  the  observation  point.  (See  Figure  1). 
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FIGURE  1.  MAGNETIC  FIELD  SOURCE  AND  OBSERVATION  POINTS 
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The  multipole  expansion  of  the  scalar  potential  is  then  obtained 
by  making  a  Taylor's  series  exapnsion  of  the  factor 


l*-M  R‘ 


(3) 


This  gives 

1/R  =  1/r  -  x;3.(l/r)  +  ~  x 'x'3^3 ^ (1/r)  +  ... 


(Aa) 


y  ^  -■-rP-  x.'  x'  ...  x.'  3  3.  ...  9.  (1/r) 

“  11  i  l2  \  li  l2  h 


( Ab) 


In  Equation  (A)  the  summation  convention  has  been  employed  wherein 
repeated  lower  case  subscripts  are  to  be  summed  from  1  to  3  with  the 
understanding  that 

X1  =  x»  x2  =  v>  x3  =  2  (5a) 

and 


3^  ■  3/3x,  3 ^  =  3/3y,  3^  =  3/3z.  (5b) 


The  summation  convention  will  be  used  throughout  this  report  unless 
otherwise  stated. 


Using  the  expansion  (Ab)  in  the  expression  for  the  scalar  poten¬ 
tial  Equation  (2)  gives 


oo 


The  integrals 


P(?') 


dV 


(6) 


(7) 


in  Equation  (6)  are  by  construction  completely  symmetric  Hth-rank  ten¬ 
sors  and  are  occasionally  referred  to  as  multipole  moments^’ 5  .  However 


(•••) 


Panofsky,  W.  K.  H.  and  Phillips,  M.,  Classical  Electricity  and 
Magnetism,  p.  15,  Addison-Wesley ,  1955. 


£ 

it  is  more  convenient  to  define  the  2  -pole  moments  in  terms  of  the  sym¬ 
metric  traceless  part  of  0.  .  Such  a  symmetric  traceless  part 

il  •••  i«. 

can  always  be  extracted  from  a  tensor  of  rank  2  or  higher  and  has  the 
following  form 


where  the  coefficients  a,b...  are  determined  by  the  equations  that 

result  from  setting  the  trace  of  Q,  ,  on  any  two  indices  equal  to 

il...  i* 

zero.  The  Greek  indices  take  on  the  values  1,2 ...l  and  the  following 
notation  has  been  used  to  denote  traces 


The  overall  factor  of  (.21)  I /2!'U  has  been  introduced  for  later 
convenience. 


(9) 


For  the  special  cases  of  t  *  1,  2,  and  3,  Equation  (8)  gives  the 
dipole  moments 


the  quadrupole  moments 


«i  ■  -f 


P (r)  Xjd  r 


<ij 


3(Qij  3  6ij  \k) 


(3xiXj 


V2> 


d3r 


(10a) 


(10b) 


and  the  octupole  moments 

^ijk  “  15^ijk  “  T^ij^mmk  +  5ik^mjm  +  6jk^imm^ 


m  j'o  (r)  [ISx^jX^  -  3('51jxk 


+  fi^x^  +  fijicx1)r^]  d^r*  (10c) 


respectively. 
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One  interesting  aspect  of  Equation  (8)  is  that  every  term  on  the 
right  side  except  the  first  involves  at  least  one  Kronecker  delta  which 
is  to  be  contracted  with  the  quantities  ...3  (1/r)  in  Equation  (6). 

These  terms  will  thus  result  in  expressions  involving  derivatives  of 
v“(l/r)  which  is  zero.  Hence,  it  follows  that 


oo 


?.=0 


s,  (1/r)  Q. 
±t  ^ 


.i 


e 


cc 


■E 


(21)1 


3.  •••3.  (1/r) 
1  Z 


(11) 


Now  by  a  process  of  induction  it  can  be  shown  that 

n,  • • -n. 


(-1)'  <20 '  *1  ‘| 
a  a  (1/r)  ,  (  » 

11  h  2  £!  T 


+  ... 


(12) 


where  the  n^,  are  the  direction  cosines  of  r;  i.e.. 


ni  =  xi^r  (13) 

and  the  missing  terms  in  Equation  (12)  all  involve  one  or  more  Kronecker 
deltas  on  the  indices  ij...i  .  Thus,  when  they  are  contracted  with  the 
traceless  multipole  moments  in  Equation  (11)  they  give  zero.  Hence, 
Equation  (11)  reduces  to 


ij>(r)  * 


(14) 


For  a  pure  2  -pole  source,  only  the  Zth  term  in  the  sum  contributes 
and  the  potential  has  the  form 


*(r) 


tf  Qi 


r 


(15) 


The  corresponding  magnetic  field  components  have  the  form 


H 


J 


(16) 
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Taking  l  derivadves  of  Equation  (16)  and  introducing  the  notation 


0  »  3  ...3  H 

J 1  *  *  * J  Jjl+1 


(17) 


results  in  the  field  derivative  equations  for  a  2  -pole  source 

D  4  =  q.  .  R.  .  . 

V,,Jl+l  1  —  — Jm-1 


(18) 


where  the  q  are  the  scaled  multipole  moments  defined  by 

„  2(1+1) 
q<  <  =  Q,  ,  /air 

1 1  * "  *  1JL  V”1*. 

is  a  (2i  +  l)-rank  tensor  defined  by 


and  R 


1 1 "  *  ’  i  i  * '  *  J 


a+i 


(19) 


d  _  2(£+l)  .  .  ji+1 

R.  .  .  ,  -  -r  S  ...d.  (n.  ...n.  /r  ) 

1r**1JlJr*,Ja+l  Ji  JA+1  ll 


(20) 


It  will  be  convenient  to  introduce  a  notation  that  obviates  the 
need  for  recording  all  of  the  subscripts.  To  accomplish  this,  paren¬ 
theses  will  be  used  to  denote  a  set  of  l  subscripted  indices  while 
square  brackets  will  be  used  to  denote  a  set  of  £+1  subscripted  indices. 
Thus, 


R(l)(j] 


£+1 


(21) 


l 

Using  this  notation,  the  2  -pole  field  Jlth-derivative  equations  take 
the  form 


°[J3  ’qU>  R(i)[j]  ’ 


(22) 


GENERAL  INVERSION  PROCEDURE 


In  this  section  a  procedure  will  be  outlined  which  is  believed  to 
be  of  general  applicability  for  the  inversion  of  Equation  (22) .  Although 
a  complete  proof  of  the  validity  of  the  procedure  in  general  has  not  been 
established,  the  procedure  is  successful  in  inverting  Equation  (22)  for 
the  special  cases  of  i  -  0,  1  and  2. 
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The  objective  is  to  use  Equation  (22)?to  determine  the  bearing 
(-nj)  and  the  scaled  moments  (q(i))  of  a  2’ -pole  source  from  a  knowledge 
of  the  field  derivatives  (D[j])  at  a  single  point.  That  this  might  be 
possible  was  first  suggested  by  W.  M.  Wynn  who  noticed  that  the  number 
of  independent  ith  derivatives  of  the  field  components  is  2S.+3,  while 
the  number  of  independent  2^-pole  moments  is  2£+l,  and  the  number  of 
independent  direction  cosines  (n^)  is  2  (except  for  a  sign  ambiguity). 
Thus,  the  number  of  independent  pieces  of  information  contained  in  the 
field  derivatives  is  precisely  the  same  as  the  number  of  independent 
unknowns.  The  number  of  independent  components  of  Dj.j  and  q(i)  follow 
from  the  fact  that  they  are  both  completely  svmmetric  and  traceless 
tensors . 

The  inversion  procedure  requires  that  two  tensors  and 

be  ^ount*  witd  tde  following  properties: 

1‘  q(i)  S(i)[j]  =  °  (23> 

2*  q(D  N(l)[j]  ^[J]  (It)  =  q(k)  (24) 


whe  re 


N(i)tj]  "  R(i)[j]  +  S(i)tj]  (25) 

'V 

N,  ....  N . . .  i  must  be  symmetric  and  traceless  in  the 
1  separate  indices  [j]  and  (k) . 

*^[ j  ]  (!)  N(i)[k]  mUSt  be  idemPotent’  1  • e •  > 

^[J](i)  N(i)UI  ^UlCm)  N(m)[k]  =  ^[j](i)  N(i)[k]  •  (26) 


By  definition  R(i)[j]  is  a  dimensionless  tensor  which  can  only  depend 
on  the  direction  cosines  n^.  Hence,  R(i)[j]  must  be  expressible  as  a 
sum  of  terms  involving  products  of  direction  cosines  and  Kronecker 
deltas.  From  Equations  (23-26)  it  then  follows  that  S^wj,  and 
N  [j](i)  must  have  similar  forms.  Thus,  one  can  write  the  most  general 
expressions  with  arbitrary  coefficients  and  then  adjust  the  coefficients 
to  satisfy  the  above  four  conditions,  thereby  obtaining  S^w.j  and 

*[J](i>  * 


Because  of  Equation  (23),  the  field  derivative  equations  can  be 
written 


D[jJ  "  q(i)  N(i)[j] 


1 


Then  contracting  on  Che  right  with  Si.  . 
Equation  (24)  ' 


yields  by  virtue  of 


q(k)  =  D[j]  ^[j](k)  ' 


(28) 


Substituting  this  back  into  Equation  (27)  results  in  the  equation 


°[i]  =  D(kj  ^[k](£)  N(i)[i]  • 


(29) 


Equation  (29)  determines  the  direction  cosines  in  terms  of  the  field 
derivatives  while  Equation  (28)  gives  the  scaled  moments  in  terms  of 
the  field  derivatives  and  the  direction  cosines.  Thus,  the  inversion 
is  complete  once  Equation  (29)  is  solved. 


Relation  (29)  involves  2£+3  equations  for  only  two  unknowns. 
Clearly  many  of  these  equations  must  be  redundant  if  there  is  to  be  a 
solution.  In  order  to  explicitly  see  the  independent  equations  in  (29) 
it  is  convenient  to  change  to  a  spherical  tensor  notationt4:>  .  For  any 
tensor  of  rank  £+1  it  is  possible  to  construct  linear  combinations  of 
the  components  which  transform  under  rotations  like  objects  of  definite 
spin  J  and  z-component  M  as  follows 


T(j) 
1 JM 


„  (j  )  JM 
K[i] 


(30) 


The  coefficients  accomplish  the  transformation  from  the  car¬ 

tesian  basis  to  the  spherical  basis  and  involve  Clebsch-Gordan  coeffi- 
cients(4,fc)  .  Since  the  construction  of  a  spin  J  object  with  z- 
component  M  from  a  set  of  £+1  spin  1  objects  (vectors)  is  usually  not 
unique,  the  superscript  set  of  £-1  indices  (j)  is  needed  to  distinguish 
the  various  ways  of  making  the  composition^)  . 


In  the  special  case  where  D[i]  is  completely  symmetric  and  trace¬ 
less  (as  it  is  in  the  current  discussion) ,  the  construction  is  unique 
so  that  Equation  (30)  is  identically  zero  unless 


J  -  J  -  £+1 
o 


Ola) 


(4 5 Rose,  M.  E. ,  Theory  of  Angular  Momentum  (Wiley,  1957). 

(6) Edmonds,  A.  R. ,  Angular  Momentum  in  Quantum  Mechanics  (Princeton 
University  Press,  1957)  Chapter  3. 

(  }Racah,  G.,  in  Ergebnisse  der  Exakten  Naturwissenschaften,  Vol.  37, 
p.  56  (Springer-Verlag,  1965). 
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and 


(j)  =  (jQ)  (31b) 

where  (j0)  is  the  unique  set  of  indices  that  gives  a  non-zero  result^'  . 
The  spherical  basis  transformation  coefficients  form  an  orthonormal 

set 


v(j>J>l  &  * 

[1]  [i]  ~  (j)(j0)  JJ’  “MM' 

and  a  complete  set 

V  jr(j)JM  _  < 

(j)JM  [k]  U]  ” 

where  the  bar  denotes  complex  conjugation  and 


(32) 


(33) 


& 


J  1,  if  (j)  -  (j ’) 

(j)  (j  ')  )  o,  otherwise 


(34) 


Equations  (32)  and  (33)  follow  from  the  completeness  and  orthonormality 
properties  of  the  Clebsch-Gordan  coefficients  themselves. 

Now  by  contracting  with  and  inserting  Equation  (33)  on  the 

right  side  Equation  (29)  becomes 


D  K(J)JM 
D[i]K[i] 


(j')J’M’ 


D<j')J’M’  -(j')J’M’  - 
lk]  K[i]  Nt«-1  (m)  N(m)[i]  *[i] 


(j)JM 

i] 

(35) 


As  pointed  out  in  Equation  (31)  this  is  a  non-trivial  result  only  when 
(j)  =  (jo)  an<i  Furthermore,  the  sum  on  the  right  side  contrib¬ 

utes  only  when  (j')  *  (j0)  and  J'  ■  J0.  Thus  the  non-trivial  part  of 
Equation  (35)  can  be  written 


T  =  T  V 
M  M’  M’M 


(36) 


^'In  the  usual  parlance  of  quantum  theory,  the  only  way  to  obtain  a 
total  spin  of  £+1  by  combining  S.+1  spin  1  objects  is  to  make  all 
the  spins  parallel. 
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or  in  matrix  notation 


T  -  TV 


where 


T  =  D  jr^Jo^o^ 
M  D[i]  K[i] 


and 


VM’M  "  KU] 


[*](»>  (rn)[ij  *[i] 


(37) 

(38) 

(39) 


In  Equation  (36)  and  subsequent  equations ,  the  summation  convention  is 
extended  to  repeated  upper  case  subscripts.  These  will  be  understood  to 
be  summed  from  1  to  2J0+1. 


The  (2J0+1)  X  (2J0+1)  matrix  Vw*M  *s  idempotent  as  can  be  seen  from 
Equations  (39),  (33),  and  (26)  and  from  the  comments  accompanying 
Equation  (31).  Thus 


v  v  =  v 
M'M"  M"M  VM'M 


(40) 


In  the  special  cases  l  ■  0,  1,  and  2  the  V-matrix  is  also  hermitian 

Vm  *  VMM’  <41> 

and  has  a  trace  given  by 


V«M  *  2J0  -  1  ■  21+1  • 


(42) 


Thus,  it  is  possible  in  these  cases  to  diagonalize  the  V-matrix  by  a 
similarity  transformation  with  a  unitary  matrix  and  obtain  diagonal 
elements  consisting  of  2fc+l  ones  and  2  zeroes.  Furthermore,  it  was 
found  in  these  special  cases  that  it  was  always  possible  to  perform  an 
appropriate  coordinate  rotation  described  by  Euler  angles  (a,6,y)  and 
thereby  diagonalize  the  associated  V-matrix.  Assuming  these  results  to 
be  valid  for  arbitrary  £,  implies  that  the  unitary  matrix  that  diago¬ 
nalizes  the  V-matrix  is  the  rotation  matrix  (°t»@»Y)  belonging  to 

the  J th  irreducible  representation  of  the  rotation  group  (Chapter  4  of 
Reference  5)(b^  . 


■b)It  should  also  be  noted  that  for  the  rotation  matrices  the  super¬ 
script  symbol  (J)  denotes  only  one  index  and  not  a  set  of  indices. 
This  and  the  two  usages  of  D  (compare  Equation  (17)  and  (45))  should 
cause  no  confusion. 
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Performing  the  similarity  transformation  on  V  in  Equation  (37) 
results  in  the  equation 


t  =  tV 


where  Vj  is  the  diagonal  form  of  V 


(43) 


(44) 


and 


Si  _  TM' 


(a,S,v)  • 


(45) 


For  the  unit  diagonal  elements  in  Equation  (43)  reduces  to  an  iden¬ 
tity.  However,  the  two  zero  diagonal  elements  give  two  non-trivial 
complex  equations  of  the  form 


S,  -  0  (46a) 

1 

Si  -  0.  <46b> 

2 


where  and  are  defined  by 
(V  )  9  (V  } 

d  ~  0  (no  sum)  .  (47) 

The  two  equations  in  Equation  (46)  are  not  independent  as  can  be 
seen  from  the  fact  that  the  usual  choice  of  phases  for  the  Clebsch- 
Gordan  coefficients  gives 


(J  ) J  M  (j  )J  -M 

v  0  0  »  /_nM  k  0  0 

K[i]  (  K[i] 


(48) 


so  that 


tm  ■  <-1>M  T-M  • 


(49) 


Furthermore, 


dm-m  (a’6’y)  -  d!m'-m  (a*s’Y) 


(50) 
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so  that 


*H  *  <-1)M  C-«  '  <51) 

Thus,  if 

tjj  =  0  (52a) 

it  follows  that 

t  M  =  0  .  (52b) 

"M1 

But  there  are  only  two  zeroes  on  the  diagonal  of  V<j  so  that  it  must  be 
concluded  that 


M 


1 


-M„ 


(53) 


and  only  one  of  the  equations  in  (46)  is  independent. 

/  g  \ 

Now  the  two  real  equations  resulting  from  the  complex  equation 


(54) 


are  the  two  independent  equations  referred  to  earlier  which  can  be 
solved  (at  least  numerically)  for  the  Euler  angles  once  the  field  deriva¬ 
tives  are  known.  After  the  Euler  angles  are  known,  the  direction 
cosines  can  be  written  immediately  (refer  to  discussion  leading  to 
Equation  (111)). 


n 

n 

n 


1 

2 

3 


-sing  cosy 

(55a) 

sing  siny 

(55b) 

cos  8 

(55c) 

(c)The  choice  of  Mi  in  Equation  (54)  seems  to  be  arbitrary.  However, 
at  least  for  lm 2,  it  is  most  expedient  to  choose  Mi  -  J0  (see 
Equation  (89)  ff.). 
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QUADRUPOLE  INVERSION 


For  a  point  quadrupole  source  (2=2)  the  potential  is  given 
Equation  (15) 

$(r)  ~  ~  Q  n^/r3  . 

Taking  the  negative  gradient  gives  the  magnetic  field 

Hk  =  I  Qij(5ninA  '  6ik  nj  '  5jk  ni)/rA  • 

Two  additional  derivatives  then  yield 
^m2k  =  ^ij  Rijk2m 


by 


with 


Dm2k  3m3A 


<ij  “  V* 

Rijk2m  =  5t63  Rijk2m  "  7(Rijk2m  +  Rijk2tn) 
+  Rijk2m  +  Rijk2iJ 


(56) 

(57) 

(58) 

(59a) 

(59b) 


(59c) 


and 


R;.,  .  =  n.n.n.  n.n 

ijk2m  i  ]t  1 


m 


R^ =  n.n.  (6.  n  +  6  n.  +  6  n  ) 
ijk2m  i  j  km  2  2m  K  k2  m 

R.  .  =  n.n„(6.  n.  +  6  n.)  +  n,  n  (6  n ,  +  6  n  ) 

ijk2m  k  2  im  j  jmi  k  m  12  j  j2i 

+  n.n  (6  ..  n.  +  6  n.) 

2  m  ik  j  Jk  i 

R^.,  „  ■  n.  (6  5  +  6  6.  +  <5.,  6  ) 

ijk2m  i  2j  km  k2  jm  jk  2m 


(60a) 

(60b) 


(60c) 


Rijkim  “  nk(5i£6jm  +  ijf6ini')  +  n£(6ikdjm  +  djk<5im) 

+  WjS  +  6jk5iP  * 


(60e) 


In  constructing  the  tensors  S(i)[j]  and  it  is  also  neces¬ 

sary  to  consider  terms  of  the  form 


R.  .  *  6  .  .n.  n.n 

ijk£m  ij  k  l  m 


Rijk8m  "  (ni6jk  "  nj6ik*  ninm  +  ^VjJl  "  nj6U)  "k11!! 
+  (ni6jm  “  nj6im)  \n£ 


(60g) 


R:..  .  *  6 . .  (n.  5 .  +  n.<5.  +  n  5.  .) 

ijkS-m  ij  k  £m  l  km  ®  ki 


(60h) 


R.  .  *  (n.6.t  -  n.6 ..  )  6.  +  (n.<5.n  -n.6..)  ■$. 

ijkdm  i  jk  j  ik  2m  i  j2  j  i SL  km 

+  (n,  <5.  -  n  6  )  fi.  .  . 

i  jm  j  im  kH 


(60i) 


All  of  the  quantities  R^,.  are  symmetric  in  (kHm)  and  either  symmetric 
or  antisymmetric  in  (ijj.  rheae  are  the  only  such  quantities  that  can 
be  constructed  as  sums  of  products  of  Kronecker  deltas  and  direction 
cosines.  Furthermore,  these  are  the  only  types  of  quantities  that  can 
occur  in  and  “[j](i)  because  of  condition  3.  Thus,  it  is  rea¬ 
sonable  to  attempt  to  construct  j  and  in  t*ie  forms 


Jijk^m 


E  •»  *“ 


ijk2m 


(61a) 


kiraij 


E  \  Ri 


ijkim 


(61b) 


After  a  lot  of  algebra  one  can  determine  the  coefficients  an  and  bn 
from  the  four  conditions  discussed  in  the  previous  section.  The  results 


a6  -  -10 


(62a) 
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and 


an  =  0,  for  n  ^  6 


(62b) 


bL  =  1/8 
b2  =  -1/1500 
b3  =  -1/24 
b4  =  1/600 
b5  =  1/60 
b6  =  1/24 

b^  =  0,  for  n  =  7,  8,  9. 
Contracting  with 


(63a) 
(63b) 
(63c) 
(63d) 
(63e) 
(63f  ) 
(63g) 


ijkim 


Rijkim  +  Sijkim 


(64) 


gives 


Nijkim  Nkimst  =  2(6is6jt  +  +  4  6ijnsnt  " 


12  V-t 


(65) 


and 


rstij  Nijkim 


C  T 

n  rstkim 


n-1 


where 

T  *  .  , .  =>  n  n  n  n,  n.n 

rstkim  r  s  t  k  i  m 

2 

T  *  n  n  n  (6,  n.  +  6.  n.  +  6,  .n  ) 

rstkim  r  s  t  km  i  Hi  k  ki  m 

T  ■  n  n  (n.  n.6^  +  n,  n  6  +  n  n  6  ) 

rstkim  r  8  k  i  tm  k  m  tH  m  t  tk 

+  n  n  (n,  n.6  +  n,  n  6  .  +  n  n„6  .  ) 

r  t  k  i  sm  k  m  si  mi  sk 

+  n.n  (n,  n„6_  +  n.  n  6  .  +  n  n.6  ’  ) 
ts  ici™  Ttmri  mi  rk' 


(66) 


(67a) 

(67b) 


(67c) 
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rstk£m  "  °rfts<6k£6tm  +  6kAi  +  Vtk* 

+  nrnt(6k£5sm  +  6km6s*  +  6ml6sk> 

+  ntn3(6k£6nn  +  Skm6r£  +  Vrk* 

TrstkJlm  =  nrInk(6s)l5tm  +  6Sm{tP  +  V5sk°tm  f  6tk6sm) 

+  nm(5sk6t*  +  6tk \l)]  +  nS[ft*'SrA»  +  6rm6t^ 

+  Wtm  +  6tk6nn)  +  Wt*  +  Vn” 

+  nt[nk('s£6rm  +  VrP  +  V6sk4rm  +  6rk6sm) 

+  n  (6  .  6  +6.6  .)] 

m  sk  ril  rk  si 

Trstk£m  “  <Vst  +  Vrt  +  Vrs>  Wm 

Trstk£m  ’  <Vst  +  Vrt  +  Vrs>  <Vkm  +  V*m  +  Vkl} 

Trstk£m  *  “kV^Vtm  +  Srt6sm  +  6st4rm* 

+  Vm(5rs6tk  +  6rt6sk  +  VW 

+  Vm(6rs6U  +  W  +  W* 


(67d) 


(67e) 

(67f) 

(67g) 


(67h) 


and 


Trstk*.m  *  ^Vktn  +  5k£6rm  +  5rk5£m^  &st 

+  (6aAm+  W  +  S*sS*J  6rt 

+  ^StkS£m  +  ^  Vkm  +  5k*.5tm^  5rs 


C1  "  C2  "  "C3  *  C6  "  4 
C4  “  C7  *  C8  "  -  12 


(671) 

(68a) 

(68b) 
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m 


(68c) 


Using  Equations  (61a),  (62a),  and  (65-68)  it  is  fairly  easy  to 
determine  that  the  four  conditions  of  the  previous  section  are  indeed 
satisfied  by  Equations  (61-63) .  As  shown  in  the  previous  section,  these 
conditions  are  sufficient  to  establish  Equations  (28)  and  (29)  which  for 
the  quadrupole  case  become 


and 


4ij 


D  Sf 
k£m  kimij 


(69) 


'■y 

Dkim  Drst  rstijNijkX,m’ 


(70) 


respectively. 

The  spherical  basis  transformation  coefficients  for  l  =  2  are 
given  by 


iki 


^  (lm1jm4  |  ij  JM)  (lm2lm3  |  llj 


m. 


m,)U  ,U  ,U  „ 
4  m^i  m^k  m^fc 


(71) 


where  (j^m ,  j  2m2 1 J  lJ  2Jm^  is  Clebsch-Gordan  coefficient^ 5  for  coup¬ 
ling  spin  ant*  spin  j2  to  obtain  spin  j  with  z-component  m;  and  U  is 
the  unitary  matrix  which  converts  the  cartesian  vector  components  to 
spherical  vector  components  (page  66  of  Reference  4). 


U  - 


5  Edmonds  notations  (Reference  5)  will  be  used  for  the  Clebsch-Gordan 
coefficients  throughout  this  report. 
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The  rows  of  U  are  labeled  from  cop  Co  boccom  by  Che  spherical  index 
m  »  1,  0,  -1  while  Che  columns  are  labeled  from  lefc  Co  righe  by  Che 
carcesian  index  i  **  1,  2,  3. 


Since  U  is  uniCary  its  inverse  is  given  by 

h  0  1  \ 

if1  -  U+  -  —  I  i  0  i  ]  (73) 

^  V  0  S5  0  / 

where  Che  rows  are  now  labeled  by  Che  carcesian  index  and  Che  columns 
by  Che  spherical  index.  FurChermore,  by  inspection  of  Equation  (72) 
it  is  easily  seen  that 


U-mi  "  ^  Umi 


(74) 


The  Clebsch-Gordan  coefficients  satisfy  an  orthogonality  relation 


^  1  2m2  M 1^ 2 ^  ^ lml^ 2m2  P 1^ 2 J 'M ' ^  (75a) 

m^m^ 

and  a  completeness  relation 


JM 

as  well  as  a  number  of  symmetry  properCies  including*5^ 

C J i~miJ 2-m2 I J 1^ 2 J-M)  “  (j1m1j2m2| jijjJH) 

,  J1+J2“J 

(J2»2J1milj2J1JM)  -  (-D  (j1n>ij2m2ljij2JM)  * 


(75b) 

(75c) 

(75d) 


(6) 


Edmonds,  A.  R. ,  Angular  Momentum  in  Quantum  Mechanics,  Princeton 
University  Press  1957,  Chapter  3. 


Using  these  relations  along  with  Equations  (72)  and  (73),  it  is  straight¬ 
forward  to  establish  a  number  of  useful  relations  for  the  spherical  basis 
transformation  coefficients. 


7jJM  '  J'M' 
ik!>  Kike 


(76a) 


jJM  - 
'iki  rs 


ir  ks  it 


(76b) 


jJ-M  .  ..J+M+l  -jJM 
Kiki  Kikl  * 


(76c) 


In  addition  it  can  be  established  that  for  the  special  case  j*2,  J=3  the 
spherical  basis  transformation  coefficients  are  completely  symmetric 
and  traceless;  i.e., 


23M  _  23M  _  23M 

iki  iik  Hcii 


(77a) 


and 


K*3"  -  0 
ikk 


(77b) 


As  discussed  in  the  previous  section  the  relations  (76b)  are  suf¬ 
ficient  to  permit  rewriting  equation  (70)  in  the  form 


T  =  T  V 
M  M'  M'M 


(78) 


where 


TM  •  "Wfi  - 


(79) 


and 


b  -23M'  jw  23M 

M'M  rst  rstij  ljkim  ^kim 


(80) 


f*’The  establishment  of  the  last  equality  in  Equations  (77a)  is  not 
quite  so  straightforward  since  it  requires  a  recombination  of  the 
three  spins  involved;  these  can  best  be  handled  by  using  Racah's  W- 
coefficients  or  Wigner's  6-j  symbols.  See  E.  P.  Wigner,  Group 
Theory  (Academic  Press,  1959),  Equation  24.24. 
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It  is  not  necessary  to  explicitly  evaluate  the  right  side  of  Equation 
(80)  in  terms  of  the  direction  cosines  although  this  could  be  done. 
Instead  it  is  only  necessary  to  establish  a  few  properties  of  the  V- 
matrix  from  the  defining  relation  (80);  namely,  that  V  has  a  trace  of  5, 
is  hermitian,  and  is  idempotent(f ^  . 

The  hermiticity  of  V  is  established  as  follows 

v  _  jf23M'  jw  „  „23M 

M'M  "  K[r]  N[r](i)  N(i)(s]  K[s] 

=  S'  C  T?  ,  r  ,  K?31?  ,  from  Equation  (66) 

Z-  n  [r]  [r][s]  [s]  M 

n-1 


EC  T?  ,  ,  ,  K?3^  ,  from  Equations  (67)  and  (77) 

„  n  [r]  [r][s]  [s] 


n=l,3,5 

Now  by  inspection  of  Equations  (67)  it  is  seen  that 


T[r][s]  =T[s][r]  *  forn“l»  3>  5 


(81) 


Thus, 


^M'M  “  £ 

n-1, 3, 5 


23M  n  — 23M'_ 

n  [s]  T[s][r]  K[r]  "  VMM’  * 


(82) 


The  idempotent  property  is  established  by  using  the  idempotent  char¬ 
acter  of  Njrw^j  [g]  along  with  Equation  (76b)  and  the  comment 

associated  with  Equations  (31) . 


(t)V  has  another  curious  property  that  follows  from  Equation  (76c). 


V 


-M-M’ 


(-1) 


M-H-T 


However,  this  relation  is  of  no  utility  in  the  present  discussion. 


r-»  —2 3M '  -.23M' 

VM'M"  VM"  M  ~  ^  K[r]  N[r](i)  N(i)[s]  [s] 


,23M 


— 2  3M "  x  V23M 

*  [t]  N[t](j)  N0)[u]  K[u] 

*  E  Vf  Via,  "cow  VT 

jJM" 

jrjJM”  a,  M  .23H 

X  K[t]  *[t](j)  N(j)[u]  K[u] 

=  K[V|  &[r](i)  N(i)[s]  ^[s]  (j)  N(j)[u]  K[u] 

=  K23M'  %  N  K23M 

K[r]  N[r](i)  N(i)[u]  K[u] 

*  VM'M  ’ 

The  trace  of  V  is  evaluated  in  a  similar  fashion 

v  =  V  K23M  Sf  N  k23m 

MM  \  [r]  N[r](i)  N(i)[s]  K[s] 


Y  rJJM  a,  23M 

^  [r]  N[r](i)  N(i)[s]  K[s] 

6[r)[s]  Vl(i)  N(i)[s]  *  N[r](i)  N(i)[r] 


=  T,  C  T?  lr  ,  . 
^  n  [r] [r] 


n=l 


(83) 


(84) 


The  quantities  T1?  . .  .  can  be  evaluated  from  Equations  (67)  with  the 
result  imrj 


V  =■  5  . 
MM 


(85) 


The  hermiticity  of  V  implies  that  V  can  be  diagonalized  by  a  simi¬ 
larity  transformation  with  a  unitary  matrix.  The  resulting  diagonal 
matrix  Vd  by  virtue  of  the  idempotent  property  of  V  can  have  only  ones 
and  zeroes  on  the  diagonal.  Finally,  the  trace  relation  (85)  requires 
that  there  be  exactly  five  ones  and  two  zeroes  on  the  diagonal  of  Vd. 
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The  position  of  the  zeroes  on  the  diagonal^can  be  established  by  con¬ 
sidering  a  special  case.  When  the  vector  r  lies  along  the  z-axis;  i.e., 

n*  -  6n  •  (86> 

The  explicit  form  of  the  V-matrix  can  be  determined  using  Equations  (67) 
and  (68)  in  the  relation 


M'M 


n=l, 3,5 


— 23M ' 
C  K 
n  rst 


^  ^3M 

rstkim  Kim 


(87) 


which  was  obtained  in  arriving  at  Equation  (82).  For  this  special  case 
one  finds  that 


V  -  Vd  =  (0,  1,  1,  1,  1,  1,  0)  (88) 

where  the  last  equality  gives  the  elements  along  the  diagonal.  From 
this  result  one  learns  not  only  the  order  of  the  elements  in  Vd,  but 
also  that  the  unitary  matrix  that  diagonalizes  the  V-matrix  corresponds 
to  a  rotation  of  coordinates  that  brings  the  z-axis  into  coincidence 
with  r.  This  implies  that  the  diagonalizing  matrix  belongs  to  a  repre¬ 
sentation  of  the  rotation  group  and  hence  is  one  of  the  rotation 
matrices  Dj^'  (cx,8,y)  where  (a,8,y)  are  the  Euler  angles  of  the  desired 
rotation.  Thus, 

[D(3) (a,6,y) VD<3) (a,6,Y>  -  Vd  =  (0,1, 1,1, 1,1,0)  (89) 

for  some  choice  of  (oi,B,y). 

(NB:  When  diagonalizing  an  hermitian  matrix  by  a  unitary  matrix, 
the  order  of  the  resulting  diagonal  elements  is  usually  immaterial. 
However,  if  the  positions  of  the  zeroes  in  Equations  (88)  and  (89)  are 
altered,  it  is  no  longer  clear  that  the  diagonalizing  matrix  is  one  of 
the  matrices  D^)(a,S,Y  )•  Thus,  to  ensure  that  the  diagonalizing 
matrix  is  one  of  these  matrices  and  not  some  other  unitary  7X7  matrix, 
it  is  expedient  to  adopt  the  order  given.) 

To  determine  the  appropriate  values  of  the  Euler  angles  the  rota¬ 
tion  is  applied  to  Equation  (78)  yielding 

t  -  tVd  (90) 


where 

t  -  TD(3)  (o.B.y) 


(91) 
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For  the  components  M  =  0,  +1,  and  +2,  Equation  (90)  is  a  trivial  iden¬ 
tity.  However,  for  M  *  +3  it  imposes  the  following  conditions  on 

U,o.>) 

t+3  -  0  .  (92) 

As  pointed  out  in  the  previous  section,  only  one  of  these  complex  equa¬ 
tions  is  independent  since  Equations  (79)  and  (76c)  imply  that 

T_m  =  <-l>M  \  •  (93) 

Thus,  it  is  sufficient  to  consider  only 

t3  -  0  .  (94) 

Using  the  notation  of  Edmonds  (Chapter  4  of  Reference  5)  for  the 
rotation  matrix,  Equation  (94)  becomes 


E  tm  ^  43)(6)el3a  =  °  •  (95) 

M 

Clearly  the  value  of  a  is  arbitrary  so  that  Equat  on  (95)  imposes  two 
real  conditions  on  the  two  angles  6  and  y. 

E  ^3)(S)  =  0  •  (96) 

M 

The  two  real  conditions  imposed  by  Equation  (96)  can  be  explicitly 
exhibited  in  a  convenient  way  by  using  Equation  (93)  and  writing 

tm  =  PMei9M  (97> 

with  and  9^  both  real.  Thus,  (suppressing  the  argument  3) 


E  v-Sff 

M  .  *.«<?  ♦  E  iv'N?* 

M*1 

■  %-of  +  Ei»M[«1<e"  +  d»>+  (-1)M  «-1(9h  * 


Taking  the  real  and  imaginary  parts  of  this  last  equality  gives  the 
desired  conditions. 


2 

<>0d03  +  £  PM  COS(0M  +  +  ("1)M  d-M3^  *  0  (99a) 

M*1 

3 

£  -  C-i)"d‘2>]  .  0  .  (99b) 

M»1 

(3) 

The  d^  (8)  are  given  in  Chapter  4  of  Reference  5. 

[ <3=m3TtWt]% 6/2  sl”3'M  2/2  •  (I00> 


Thus, 


d(3) 

33 

(6)  *  cos6  6/2 

(101a) 

d<3) 

d23 

(6)  =  -  v^T  cos3  6/2  sin  6/2 

(101b) 

d(3) 

d13 

(6)  ■  v/l5  cos^  6/2  sin2  6/2 

(101c) 

d(3) 

d03 

(6)  -  -2/5  cos3  6/2  sin3  6/2 

(lOld) 

d(3) 

d-13 

(6)  *  /15  cos2  6/2  sin^  6/2 

(lOle) 

d(3) 

-23 

(6)  ■  -/6~ cos  6/2  sin3  6/2 

(lOlf) 

d(3) 

-33 

(6)  -  sin6  6/2  . 

(101g) 

Defining 

djJ3)  (8)  +  (-1)M  d^3j  (6) 

(102a) 

and 

\  ■  d<2)  t»  -  <-»M  d<2>  (9). 


(102b) 


one  finds 


*0  '  ’  -  4  Sl”3  6 

(103a) 

/l5  3 

A1  =  r—  sin  S  cot  0 

(103b) 

A2  =  ~  sin3  0(1  +  2  cot2  0) 

(103c) 

A3  =  ^  sin3  S  cot  0(3  +  4  cot2  0) 

(103d) 

Bx  =  ~  sin2  0 

(103e) 

6  2 

B2  =  — ^  sin  0  cot  0 

(103f ) 

B3  =  ^  sin2  0(1  +  4  cot2  0)  . 

(103g) 

Using  Equations  (102-103)  and  defining 

X  -  cot  0, 

Equations  (99),  after  some  simplification,  become 
sin3  0[p3x(3  +  4x2)  cos(03  +  3y) 

-  P2(l  +  2x2)  cos(92  +  2y) 

+  /~15  pxx  cos (0^  +  y)  -/$  Pq]  -  0 

sin20  Ip3(1  +  4x2)  sin(63  +  3y) 

-2/6*  P2x  sin(©2  +  2y) 

+  1^15  p^  sin(0^  +  y)]  •  0  .  (105b) 

Solution  of  these  two  equations  for  g(or  x)  and  y  then  permits  the  com¬ 
pletion  of  the  inversion  process.  Unfortunately  these  equations  are 
transcendental  and  do  not  seem  to  permit  a  simple  solution.  They  are 


(104) 


(105a) 
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probably  best  solved  numerically  once  the  values  of  and  are  known. 
Only  solutions  in  the  ranges 

0  <  0  <  n  (106a) 

0  <_  y  <  2tt  (106b) 

need  be  determined.  Once  the  solutions  of  Equations  (105)  are  known  the 
direction  cosines  nj  can  be  easily  determined. 

The  spherical  basis  components  Ym  of  n.i  are  given  according  to 
Equations  (30)  and  (71)  by 


Y 

m 


U  ,n. 
mi  i 


(107) 


In  the  special  frame  where  n^^  *  6^3  (denoted  by  a  prime)  this  gives 


Bv  Equation  (91)  Y'  and  Y 
mm 


Y'  *  U  _  =  <5 
m  m3  mo 

are  related  by 


(108) 


Ym  *  Ym»  Dm'm  (a .V)  •  (109) 

m  m  mm 

Solving  Equations  (107-109)  for  n^  gives 

n.  -  uTilD^^.e.Y)]"1  .  (no) 

1  im  om 

Explicit  evaluation  of  Equation  (110)  using  Edmond's  expression  for 
d(1)  (c*,B,y)  gives 


n^  ■  -sin  B  cos  y 
n2  -  sin  6  sin  y 
n^  ■  cos  0  . 


(111a) 

(111b) 

(111c) 


It  i3  now  straightforward  to  evaluate  the  scaled  moments  q^j  from 
Equation  (69),  thus  completing  the  inversion  process  for  the  quadrupole 
source. 

2  4 

NB:  It  should  be  noted  that  the  terms  involving  Rijkim  and  Rijkim 
do  not  contribute  to  the  sum  in  Equation  (69)  because  of  the  traceless 
character  of  Furthermore,  the  evaluation  of  the  remaining  sums 

can  be  greatly  simplified  by  taking  advantage  of  the  symmetry  of  D.  .  . 


26 


t 


SUMMARY  OF  THE  QUADRUPOLE  INVERSION  RECIPE 


Once  the  values  of  the  independent  second  derivatives  (for  example 

^111>  ^112>  ^113>  ^122»  d123»  0222,  ®223^  are  known,  the  spherical  com¬ 
ponents  can  be  computed  from  Equation  (79).  After  considerable  simpli¬ 


fication  these  give 

T3  =  -/2/2[Dm  -3D122  +  i(3D112  -  D^)]  (U2a) 

T2  =  /3/2(D113  -  D223  +  2i  D123]  (112b) 

Tj_  =  /30/4[Dm  +  D122  +  i(D112  +  D222)]  (112c) 

To  =  -^0/2[Dn3  +  D223]  .  (112d) 

The  moduli  pM  and  arguments  9^  can  then  be  computed  from  the  real  and 
imaginary  parts  of  TM  in  the  usual  way. 

v 

PM  =  [(Re  Tm)2  +  (Im  TM)2]2  (113a) 

sin  eM  -  Im  Tm/pm,  cos  8m  *  Re  Tm/pm  .  (113b) 


These  can  then  be  substituted  into  Equations  (105)  and  the  resulting 
equations  solved  for  6  and  y  in  the  intervals  specified  in  relations 
(106). 

For  convenience,  the  pertinent  equations  are  repeated  here. 
sin38[p3x(3  +4y2)  cos(e3  +  3y) 

-t/6~ p2(l  +  2x2)  cos(e2  +  2y) 

+  /l5  p^xcosO^  +  y)  -  i/5~ pQl  (105a) 

sin28[p3(l  +  Ax2)  sin(93  +  3y) 

-2  «/6~  p2x  sin(e2  +  2y) 

+  t/15  p3  sln^  +  y)] 


(105b) 


0  ±  6  £  n  (106a) 

0  <  y  <  2tt  .  (106b) 

Finally,  the  direction  cosines  are  given  by  Equations  (111) 

n^  =  -sin  S  cos  y  (111a) 

n2  *  sin  0  sin  y  (111b) 

n^  ■  cos  6  (111c) 

and  the  scaled  quadrupole  moments  are  given  by  Equations  (69),  (61b), 
(60)  and  (63).  After  much  simplification  the  latter  give(ti 


qij  "  10  nmDmij  "  B^iVAti  +  “iVj^kiJ* 


+  4(n.n,  +  4  6.  .)  n.  n.n  0,  „ 

8  i  j  3  ij  k  £  ro  k£m 


(114) 


For  purposes  of  checking  the  results,  the  field  derivatives  can  be 
regenerated  using  Equations  (58)  through  (60).  When  combined  these  can 
be  simplified  to  give 

°k£m  “  5{7^9nkn£nm  '  (6kmn£  +  6 lm\  +  5klntn^  ninjqij 

+  2[(-7nknt  +  6^)  n^  +  (- 7 +  *km) 


+  (~7n£nm  +  V 

+  2(nkq£m  +  n£qkm  +  Vk«)}  * 


(115) 


U) 


See  comments  following  Equation  (111) . 


APPENDIX  A 


DIPOLE  INVERSION  RECIPE 
USING  THE  SPHERICAL  TENSOR  PROCEDURE 


From  the  independent  field  derivatives  (D,.,  D10,  D..,,  D,0,  D.J, 


evaluate  the  spherical  components  T^. 


T2  *  2(D11  “  °22  +  2lD12) 


T1  *  -(D13  +  iD23> 

To  *  <D  +  D22> 


(Ala) 

(Alb) 

(Ale) 


The  moduli  and  arguments  9^  can  then  be  evaluated  in  the  usual  way 


h 


PM  =  [(Re  V2  +  (Im  V2] 

sin  0M  =  Im  cos  0M  »  Re  TM/pM  . 


(A2a) 


(A2b) 


The  angles  6  and  y  are  then  obtained  by  solving  the  following  system  of 
equations 


sin2S[2p2(l  +  2y2)  cos(02  +  2y) 

-4p^xcot(8^  +  y)  +  /fe”  p  ]  =  0  (A3a) 

sin  8[p2x  sin(92  +  2y)  -  p^  sin(91  +  y)]  *  0  (A3b) 

with 

0  <  s  <  TV  (A4a) 

0  <_  y  <  2it  .  (A4b) 

The  direction  cosines  are  given  by 

n^  »  -sin  g  cos  y  (A5a) 

n2  ■  sin  6  sin  y  (A5b) 

n^  ■  cos  8  .  (A5c) 


The  scaled  dipole  moments  are  given  by 

’t  ■  I  \Du  - 1  VAV  • 

For  checking  purposes  the  field  derivatives  are  regenerated 
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